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We  study  the  thermoelectric  transport  properties  of  a  single  electron  transistor.  We  describe  a  two-level 
quantum  dot,  connected  to  right  and  left  leads,  employing  the  single  impurity  Anderson  model  with  local 
finite  electronic  correlation.  Using  the  linear  response  theory,  we  compute  the  thermoelectric  transport 
coefficients.  We  calculate  the  thermoelectric  properties  employing  the  Green's  functions  calculated 
within  the  atomic  approach,  but  for  simplicity,  we  only  consider  the  electronic  contribution  and  neglect 
the  phononic  contribution  to  the  thermal  conductance.  In  the  single  electron  transistor  this  is  not  so 
drastic,  because  the  phononic  contribution  can  be  minimized,  isolating  the  quantum  dot  from  the 
electrodes  with  two  tunneling  barriers,  whose  material  can  be  appropriately  chosen  to  produce  a  low 
phononic  contribution.  We  show  that  the  best  dimensionless  thermoelectric  figure  of  merit  for  the  single 
electron  transistor,  occurs  in  the  weak  coupling  regime,  at  temperatures  well  above  the  Kondo 
temperature. 

©  2014  Elsevier  Masson  SAS.  All  rights  reserved. 


1.  Introduction 

The  development  of  semiconductor  technology  in  the  1950's, 
with  the  discovery  that  heavily  doped  semiconductors  and  the 
compound  B^Tes  are  good  thermoelectric  materials,  gave  rise  to 
the  modern  thermoelectric  industry,  with  the  development  of  new 
commercial  power  generation  and  cooling  systems.  During  the  next 
three  decades,  with  little  interest  from  the  scientific  community, 
the  thermoelectric  field  developed  slowly.  But  at  the  beginning  of 
the  1990's,  after  the  publication  of  two  important  papers  by  L.D. 
Hicks  and  M.  S.  Dresselhaus  [1],  the  situation  changed  completely. 
They  suggested  that  the  quantum  confinement  present  in  low 
dimensional  systems  could  be  employed  for  engineering  the  ther¬ 
moelectric  properties  of  those  materials.  In  the  wake  of  these  pa¬ 
pers,  an  intense  effort  was  made  to  explore  the  new  routes  opened 
by  the  experimental  achievements  in  nanoscopic  systems  research, 
such  as  quantum  dots,  nanowires,  quantum  wells,  superlattices, 
and  more  recently  nanocomposite  materials.  A  critical  discussion  of 
the  development  of  this  important  area  can  be  found  in  some 
recent  reviews  [2.3, 4, 5, 6]. 
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During  the  last  10  years,  several  experimental  setups  for 
exploring  the  thermoelectric  properties  of  nanoscopic  systems 
[7-11]  have  opened  up  the  possibility  of  using  such  systems  for 
controlling  the  energy  transport  on  a  microscopic  scale.  In  a  recent 
experiment,  a  gate-controlled  single-electron  refrigerator  was  ob¬ 
tained,  in  which  the  Coulomb  correlation  controls  the  heat  transfer 
and,  the  refrigeration  process  in  a  mesoscopic  system  [12].  In 
another  interesting  paper,  Popescu  et  al.  [13],  investigated  the 
fundamental  dimensional  limits  for  thermodynamic  machines  and 
showed  that  it  was  possible  to  construct  self-contained  re¬ 
frigerators  (i.e.,  not  requiring  external  sources  of  work)  consisting 
of  only  a  small  number  of  qubits  and/or  qutrits.  Another  paper 
related  to  our  calculations  was  the  study  of  thermoelectric  effects  in 
molecular  junctions  by  Pramod  Reddy  et  al.  [14],  who  trapped 
molecules  between  two  gold  electrodes  with  a  temperature  dif¬ 
ference  across  them.  From  the  theoretical  point  of  view  we  can 
refer  two  recent  papers  employing  the  numerical  renormalization 
group  [15,16]. 

A  strong  motivation  for  studying  the  single  electron  transistor 
(SET)  is  that  this  system  is  the  experimental  realization  of  the  single 
impurity  Anderson  model  (SIAM)  for  finite  electronic  correlation  U. 
The  SIAM  was  experimentally  realized  by  the  Goldhaber-Gordon 
group  [17]  with  complete  control  over  all  the  parameters  of  the 
model.  They  measured  the  electric  conductance  of  a  SET  and 
showed  its  universal  character. 
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Nomenclature 

r 

energy  unit,  J 

A 

Anderson  parameter,  J 

nF 

Fermi-Dirac  distribution, 

Constants 

I<e 

electronic  thermal  conductance,  W  nrr2  K-1 

kB 

Boltzmann  constant  =  1.38  x  10  23  m2  Kg  s  2 

Kph 

phonon  thermal  conductance,  W  m-2  I<-1 

e 

electron  charge  =  1.6  x  10-19  C 

Wi,/2 

transport  coefficients, 

Lo 

Lorenz  number,  W  Q  I<-2 

Symbols 

S 

Thermopower,  V  K-1 

e 

electron  charge,  C 

T 

temperature,  I< 

chemical  potential,  J 

ZT 

dimensionless  thermoelectric  figure  of  merit, 

V 

hybridization,  J 

G 

electrical  conductance,  C2  Kg-1  itT2  s'1 

U 

correlation  energy,  J 

1 0 

energy,  J 

Eqp 

quantum  dot  energy,  J 

Q 

thermodynamical  potential,  J 

D 

semi-width  conduction  band,  J 

The  dimensionless  thermoelectric  figure  of  merit  ZT  indicates 
the  system  performance.  It  constitutes  a  measure  of  the  usefulness 
of  materials  or  devices  to  be  employed  for  thermopower  generators 
or  cooling  systems,  and  is  defined  by. 

zt  =  s2tg/k,  (1) 

where  G  is  the  electrical  conductance,  S  is  the  thermoelectric  po¬ 
wer,  I<  is  the  thermal  conductance  and  T  is  the  absolute  tempera¬ 
ture.  If  ZT  >  1,  the  system  is  considered  to  be  a  possible  candidate 
for  technological  applications.  In  ordinary  metals,  the  thermal  and 
electrical  conductances  are  related  to  the  same  scattering  pro¬ 
cesses,  with  only  weak  energy  dependence.  This  is  the  main  reason 
why  metals  show  lower  dimensionless  thermoelectric  figure  of 
merit  ZT  values.  After  the  discovery  in  the  fifties,  that  the  alloys  of 
Bi2Te3  exhibit  ZT  -  1,  near  room  temperature,  this  compound  has 
dominated  the  whole  field  of  thermoelectric  materials  up  until 
today.  All  efforts  to  increase  the  ZT  values  have  resulted  in  few 
advances,  due  to  the  interdependence  between  the  physical  prop¬ 
erties  associated  with  ZT.  A  summary  of  the  evolution  of  ZT  values 
can  be  found  in  Fig.  1  of  the  C.  J.  Vineis  et  al.  review  [4].  However, 
nanoscopic  systems  open  up  new  possibilities  for  increasing  ZT, 


Fig.  1.  (a)  Schematic  picture  of  the  SET.  A  quantum  dot  embedded  into  conduction 
leads  and  (b)  Details  of  the  eigen  states  structure  of  the  atomic  solution  generated 
through  the  atomic  approach  in  the  Kondo  regime. 


mainly  due  to  the  level  quantization  and  the  Coulomb  interaction, 
which  could  lead  to  important  changes  in  the  thermoelectric 
properties  of  the  system. 

In  the  present  paper,  we  show  that  the  violation  of  the 
Wiedemann-Franz  law  (WF)  within  the  temperature  interval 
kBT  =  [1  — 5]T,  where  kB  is  the  Boltzmann  constant,  is  associated 
with  a  strong  increase  in  the  dimensionless  thermoelectric  figure  of 
merit  ZT.  T  is  an  energy  unit  associated  with  the  energy  width  of 
the  conduction  band.  We  assume  a  flat  density  of  states  of 
magnitude  pc  =  1/2 D  per  spin  channel  for  both  leads,  where  we 
employ  the  energy  unit  T  =  1  with  D  =  100r  being  the  half 
bandwidth  of  each  lead.  To  simplify  the  notation  we  set.  kB  =  1. 

The  conductance  unitary  limit  is  experimentally  attained  in  real 
quantum  dot  (QD)  systems,  at  temperatures  on  the  order  of  10  mK. 
W.  G.  van  der  Wiel  et  al.  [18]  obtained  the  unitary  limit  in  a  semi¬ 
conductor  QD  at  T  =  15  ml<  with  an  applied  magnetic  field  of 
B  =  0.4T,  and  A.  V.  Kretinin  et  al  19]  obtained  the  same  limit  in  a 
suspended  nanowire-based  QD  device,  at  temperature  T  =  10  mK. 
In  their  setup,  there  is  no  applied  magnetic  field.  Since  in  our  units, 
the  unitary  limit  of  the  conductance  was  attained  at  T  =  lCT3r 
(these  results  are  not  shown  here),  we  estimate  r  -  10  K,  and  since 
we  are  studying  the  temperature  region  T  =  [  1  — 5]T  where  charge 
fluctuations  becomes  important,  we  estimate  that  in  real  systems, 
those  temperatures  will  occur  approximately  within  the  range  of 
T  —  [10-50]  K,  which  is  easily  attained  in  a  liquid  helium  cryostat. 

In  a  general  way,  G.D.  Mahan  and  J.  O.  Sofo  20]  discussed,  the 
conditions  that  a  device  must  fulfill  in  order  to  produce  the  best 
thermoelectric  efficiency.  They  investigated  what  type  of  electronic 
structure  provides  the  largest  dimensionless  thermoelectric  figure 
of  merit  for  thermoelectric  materials.  They  determined  that  a  nar¬ 
row  energy  distribution  of  the  carriers  was  needed  in  order  to 
produce  a  large  value  of  ZT.  In  the  present  paper,  we  show  that  this 
general  condition  is  fulfilled  by  the  SET  at  the  weak  coupling  limit 
(low  hybridizations  values)  and  within  the  temperature  range 
T  —  [10-50]  K,  above  the  Kondo  temperature.  The  regions  where 
ZT  exhibits  high  values  corresponds  to  the  transition  from  the 
empty  dot  to  the  intermediate  valence  and  the  transition  from  the 
singly  to  the  doubly  occupied  states.  We  show  that  this  is  the  most 
interesting  region,  because  ZT  attains  strong  values. 

Considering  the  electronic  and  phononic  contributions  sepa¬ 
rately,  ZT,  can  be  written  as. 

ZT  =  S2TG/  (l<e  +  Kph) .  (2) 

where  I<e  and  KPh  are  the  electronic  and  the  phononic  contribution 
to  the  thermal  conductance,  respectively. 
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In  low  dimensional  structures,  such  as  films,  wires,  and  super¬ 
lattices,  the  large  values  of  the  phononic  contribution  attained  by 
JCph  is  responsible  for  low  ZT  values  [21].  In  the  SET,  the  phononic 
contribution  to  thermal  conductance  can  be  minimized,  because 
the  QD  is  isolated  from  the  electrodes  by  two  tunneling  barriers, 
whose  material  can  be  appropriately  chosen  to  produce  low  Kph 
values  or  by  alloying  the  barriers  material  in  order  to  introduce 
atomic  scale  defects  and  scatter  phonons  [4].  However,  since  we  are 
working  at  temperatures  within  the  interval  T=  [1  -5]r,  we  believe 
that  the  increase  of  the  phononic  contribution  is  important  for 
diminishing  the  strong  ZT  increase  due  to  electrons. 

Taking  into  account  that  usually  a  SET  is  basically  a  lithographic 
nanostructure  obtained  over  a  GaAs(i_X)Al*-GaAs  quantum  well 
[7-9,17  ,  where  electrons  and  holes  are  confined  in  a  two  dimen¬ 
sional  gas,  and  the  substrate  that  “supports”  the  two  dimensional 
gas  consists  of  two  GaAs(i_X)Alx  bulk  pieces,  each  one  with  a  Debye 
temperature  To  —  300  K  22],  and  since  our  estimation  of  the 
temperature  interval  is  T  =  [1  -5]T  -  [10-50]  K,  we  confirmed  that 
the  temperature  values  are  below  the  value  of  the  Debye  temper¬ 
ature  of  the  substrates,  as  a  consequence,  we  expected  that  the 
phononic  contribution  of  the  substrates  would  be  low.  Recently, 
new  research  has  been  published  that  takes  into  consideration  the 
phononic  contribution  to  the  thermal  conductivity.  That  most 
similar  to  the  present  paper  is  that  of  Jie  Ren  et  al.  [23].  However,  in 
order  to  include  phonons  we  need  to  reformulate  the  atomic 
approach  which  can  make  the  calculations  very  complicated,  and 
for  the  sake  of  simplicity,  we  do  not  consider  the  phononic 
contribution  to  the  thermal  conductivity  Kp h  in  the  present  paper, 
but  rather  we  indicate  its  main  effect  on  the  thermoelectric 
properties. 

We  organize  this  paper  in  the  following  way:  In  Section  2,  we 
define  the  model  in  terms  of  the  X  Hubbard  operators.  In  Section  3, 
we  present  a  brief  review  of  the  main  equations  employed  in  the 
calculation  of  the  Green's  functions  by  the  atomic  approach.  In 
Section.  4,  we  discuss  the  thermoelectric  properties.  In  Section  5, 
we  present  the  results  of  our  research  and  discuss  their  physical 
consequences.  Finally,  in  Section  6,  we  present  a  summary  of  the 
main  results  obtained  as  well  as  the  conclusions  of  the  paper. 


2.  Model  and  theory 


In  Fig.  la,  we  present  a  pictorial  view  of  the  SET  considered  in 
the  present  paper.  In  Fig.  lb,  we  represent  the  eigen  states  of  the 
atomic  solution,  which  are  the  “seed”  employed  to  calculate  the 
Green's  functions  of  the  atomic  approach.  Those  Green's  functions 
have  the  general  form  [24]: 


i=\  °J  ~  Ul 


(3) 


where  Cl  is  the  thermodynamic  potential,  the  poles  u*  in  the  de¬ 
nominator,  are  the  energy  eigenvalues  associated  with  the  energy 
levels,  and  the  residues  m,  in  the  numerator,  are  the  occupation 
numbers  associated  with  the  atomic  transitions.  In  the  diagram,  we 
indicate,  transitions  with  null  residues  ( m*  =  0),  with  dashed  lines 
and  transitions  with  non  null  residues  (m,  #  o)  with  solid  lines,  but 
as  the  temperature  or  the  other  parameters  of  the  model  change, 
the  residues  of  the  Green's  functions  change  as  well  and  the  dotted 
lines  can  become  solid  lines  and  vice-versa.  The  resultant  eigen 
state  structure  of  the  QD  corresponds  to  the  atomic  solution  of  the 
periodic  Anderson  model  (PAM),  (For  more  details  cf.  Fig.  20  and 
the  appendix  of  Ref.  [24]). 

This  particular  eigen  state  structure  corresponds  to  the  Kondo 
limit  at  a  temperature  T,  below  the  Kondo  temperature  TK.  The 


Kondo  singlet  ground  state  is  given  by  the  level  11,  and  the  tran¬ 
sitions  represented  by  the  full  lines  produce  non  null  residues  for 
the  atomic  Green's  functions,  which  give  rise  to  the  Kondo  peak. 
The  transitions  represented  by  dashed  lines  vanish  in  the  Kondo 
region,  but  become  important  in  the  temperature  region,  where  we 
obtain  high  ZT  values.  This  QD  eigen  state  structure  is  not  static,  but 
rather  changes  with  temperature  and  the  QD  localized  energy  level 
EfCT,  follows  the  different  regimes  of  the  SIAM:  empty  dot,  inter¬ 
mediate  valence,  Kondo,  and  double  occupation. 

The  Hamiltonian  of  the  system  can  be  written  as. 

w  =  E  E  vlAa*  +  E  M-  +  (£a  +  u)*dd] 

a-L,R  k,a  "  (4) 

+  E  EV«(*U,.+<«X*). 

a=L,R  k,a 

where  the  first  term  represents  the  leads,  characterized  by  free 
conduction  electrons  (c-electrons)  to  the  right  (R)  and  to  the  left  (L) 
of  the  QD  (Fig.  la).  The  second  term  describes  the  QD  defined  by  a 
two-level  structure:  one  localized  (/-electrons)  bare  level  EfG  =  EqD 
and  the  local  Coulomb  interaction  of  the  QD,  characterized  by  the 
double  occupation  level  at  Eqd  +  hi  [24,25].  These  localized  levels 
are  expressed  in  the  representation  of  Hubbard  operators,  which 
are  convenient  for  working  with  correlated  local  states  and  are 
defined  in  general  by  Xp>ab  =  |p,axp,jb|.We  would  like  to  mention 
that  the  simple  two-level  structure  of  the  isolated  QD,  is  changed 
when  this  QD  is  allowed  to  interact  with  the  conduction  electron 
band.  This  simple  structure  is  modified  by  the  interaction  and  the 
atomic  approach  changes  it,  into  the  “atomic  eigen  state  structure” 
represented  in  Fig.  lb.  The  third  term  corresponds  to  the  tunneling 
between  the  embedded  dot  and  the  left  (I)  and  right  (R)  semi¬ 
infinite  leads.  The  amplitude  Va  is  responsible  for  the  tunneling 
between  the  QD  and  the  lead  a.  For  simplicity,  we  assume  sym¬ 
metric  junctions  (i.e.  Va  =  Vi  =  Vr  =  V)  and  identical  leads  con¬ 
necting  the  QD  to  the  quantum  wire. 

We  solve  the  Hamiltonian  employing  the  atomic  approach  for 
the  finite  correlation  U  case.  This  method  was  developed  in  earlier 
papers  [24,25],  and  here  we  present  a  short  review  of  it  in  Section  4. 
The  atomic  approach  is  able  to  capture  the  low  temperature 
physics,  dominated  by  the  Kondo  effect  as  well  as  the  region  T  >  T, 
where  charge  fluctuation  processes  are  dominant.  In  that  temper¬ 
ature  region,  the  dimensionless  thermoelectric  figure  of  merit  ZT 
attains  values  greater  than  the  unit,  in  the  transition  from  the 
empty  dot  to  the  intermediate  valence  regions. 

3.  The  atomic  approach  for  the  single  impurity  Anderson 
model  (SIAM) 

In  previous  papers  [24,25],  one  of  us  developed  a  simple 
method,  that  employs  the  atomic  solution  of  the  periodic  Anderson 
model  as  a  “seed”,  for  generating  approximate  Green's  functions  for 
describing  the  Kondo  physics  of  strongly  correlated  electron 
systems. 

The  exact  GF  for  the /electrons,  which  is  valid  both  for  the  PAM 
and  the  SIAM,  can  be  written  in  the  form  of  a  Dyson  equation  [24] 

=  Mtf-  (I  —  A^)-1,  (5) 

where  Ac  =  W^M^  and  Ma  are  the  exact  effective  cumulants. 
Inverting  Eq.  (5),  we  obtain 

M„  =  (i  +  g£-W„)-1.g£,  (6) 

and  for  an  impurity  located  at  the  origin  Wa  is  given  by. 
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WT(z)  =  |V|2n(z)I, 

(7) 

wi(z)  =  |V|2<p|(z)I', 

(8) 

where 

-0  0  •  “A 

(9) 

For  a  rectangular  band  with  half-width  D  in  the  interval  [-D,D], 
we  have. 


Hamiltonian  (For  more  details  cf.  the  appendix  of  our  reference 
[24]).  From  this  equation  we  obtain  the  exact  atomic  cumulant. 

Ma; 

Mf=  (i  +  g^-W (17) 
where 

W°T  (z)  =  |A|2«p°  (Z)  ;  W°  (Z)  =  |A|2«p°  (z)I,  (18) 

and. 


<p4z)  =  2£jln 


z  +  D  -F  /x\ 
z-D  +  n) 


(10) 


(19) 


where  /jl  is  the  chemical  potential. 

Since  the  calculation  of  the  exact  effective  cumulant  is 
equivalent  to  obtaining  the  exact  Green's  functions  [24,25],  we 
introduce  an  approximation  that  consists  of  substituting  the  exact 
effective  cumulant  in  the  Green's  functions  with  the  approximated 
ones,  obtained  from  the  atomic  solution  of  the  model,  which  is 
associated  with  the  transitions  described  by  the  Hubbard  operators. 
We  employ  the  index  Ix  =  1,2, 3, 4,  defined  in  Table  1,  to  characterize 
these  X  operators: 

In  the  atomic  approach,  those  approximated  GF's  become. 

’  i-iviVms13 

with  R 13  =  mnm33-mi3m3i  and.  S13  =  *7711+77733+77713+77731 


This  equation  corresponds  to  the  zero-width  band  located  at  eo, 
namely  the  bare  conduction  Green's  function.  The  atomic  cumulant 
Mf,  (Eq.  (17))  is  obtained  by  replacing  Eqs.  (18)-19)  in  Eq.  (17).  The 
full  Green's  functions  g£  are  calculated  by  substituting  the  cumu- 
lants  obtained  in  this  way  in  Eqs.  (12)  and  (13).  This  procedure 
overestimates  the  contribution  of  the  c  electrons,  because  we 
concentrate  them  at  a  single  energy  level  e0 ,  and  to  moderate  this 
effect  we  replace  V2  in  Eq.  (18),  with  the  Anderson  parameter 
A  =  7tV2pc ill)  =  7rV2/2D. 

We  should  stress  that  in  the  low  temperature  regime  the  atomic 
approach  imposes  the  fulfillment  of  the  Friedel  sum  rule,  but  for 
temperatures  above  the  Kondo  temperature  T/<,  as  happens  in  the 
present  paper  for  T  -  F,  we  do  not  apply  this  procedure;  we  only 
put  the  conduction  atomic  level  of  Eq.  (19)  at  the  chemical  potential 
position,  eo  =  /jl. 


d (i3) 

with  R24  =  77122*7144  -  7712477142  and  S24  =  71122  +  *7144  +  *7124  +  *7142  and. 


are  the  atomic  cumulants  calculated  from  the  atomic  Green's 
functions  (For  more  details  cf.  section  III  and  the  appendix  of  our 
Ref.  [24]).  In  the  same  way,  we  can  obtain  the  conduction 
G^c(k,k',iw)  and  the  cross  GCJ (k,  ioj)  Green's  functions. 

The  exact  atomic  Green's  function  g£,at(z)  of  the  atomic  problem 
satisfies  a  Dyson  equation  of  the  same  form  as  Eq.  (5),  but  now  in 
terms  of  the  atomic  cumulants  given  by  Eqs.  (14)  and  (15). 

g+  =  Mat(I-W°M If)'1.  (16) 

It  should  be  stressed  that  g£,at  is  known  and  it  is  obtained  from 
the  exact  analytical  diagonalization  of  a  16  x  16  matrix,  which 
represents  the  atomic  solution  of  the  periodic  Anderson 


Table  1 

Representation  of  the  possible  transitions  present  in  the  finite  U  atomic  SIAM 
Hamiltonian.  Ix  =  1,3  destroys  one  electron  with  spin  up,  and  Ix  =  2,4  destroys  one 
electron  with  spin  down. 


(ii) 


Ix 

1 

2 

3 

4 

a  =  (A,  a) 

(0.T) 

(0,1) 

(T,<0 

4.  Thermoelectric  properties 


For  ordinary  metals,  the  Wiedemann-Franz  (WF)  law  states  that 
the  relation  between  the  electronic  contribution  to  the  thermal 
conductance  I<e  and  the  product  of  the  temperature  T  and  the 
electrical  conductance  G. 


WF  = 


Keen 

TG(T)  ’ 


(20) 


is  independent  of  temperature  and  assumes  a  universal  value  given 
by  the  Lorenz  number  [16],  L0  =  (7r2/3)(/<B/e)2,  where  I<b  is  the 
Boltzmann  constant  and  e  is  the  electron  charge.  In  our  calcula¬ 
tions,  we  employ  WF  defined  in  terms  of  L0  units. 

Combining  Eqs.  (1)  and  (20)  we  can  write  ZT  as. 

ZT  =  S2TG  j  Ke  =  (21) 

This  equation  shows  that  we  can  lead  to  cause  an  effective 
growth  of  ZT  in  regions  where  WF  <  1  and  at  the  same  time  S  >  1. 
We  found  this  optimum  combination,  within  the  temperature  in¬ 
terval  T  =  [  1  — 5]T,  at  the  weak  coupling  limit  of  the  model,  when 
the  hybridization  is  low.  In  those  regions,  we  obtain  higher  ZT 
values. 

The  thermoelectric  properties  are  calculated  following  the 
standard  textbook  definitions  [26,27].  In  this  way,  the  electrical 
conductance  G(T),  the  thermal  conductance  Ke(T )  and  the  ther¬ 
mopower  (Seebeck  effect)  are  defined  by  the  relations 
respectively. 


C(T)  =  e2I0(T), 


(22) 
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Fig.  2.  Thermopower  S  vs.  the  gate  voltage  £qd  for  different  hybridizations  V  values,  at 
temperature  T  =  T. 


Fig.  4.  Dimensionless  thermoelectric  figure  of  merit  ZT  vs.  gate  voltage  %D  for 
different  hybridizations  V  and  temperature  T  =  T.  In  the  inset  a)  of  the  figure  we  plot 
the  parameter  In  the  inset  b)  of  the  figure,  we  plot  for  V  =  4.0r;  the  ZT  values,  the 
Wiedemann-Franz  law  WF  and  the  thermopower  S. 


Ke(T) 


and. 


h(T) 

I0(T)  * 


(23) 


(24) 


To  calculate  the  transport  coefficients  /0(T),  7i(T),  and  72(T),  we 
follow  the  paper  by  Dong  and  X.  L.  Lei  28],  which  derived  the 
particle  current  and  thermal  flux  formulas,  through  an  interacting 
QD  connected  to  the  leads,  within  the  framework  of  the  Keldysh 
non-equilibrium  Green's  functions  (GF).  The  electric  and  thermo¬ 
electric  transport  coefficients  were  obtained  in  the  presence  of  the 
chemical  potential  and  temperature  gradients  with  the  Onsager 
relation  in  the  linear  regime  automatically  satisfied.  The  transport 
coefficients  consistent  with  the  thermoelectric  general  formulas 
derived  earlier  are  given  by. 


dnF(u,  T)\ 
dco  ) 


WnT(w)dw, 


(25) 


where  nF(e,T )  =  1/(1  +e^_/x)/kBT)  is  the  Fermi-Dirac  distribution, 

ra(oj)  =  (V2/A)2|g0oH)|  is  the  transmittance  for  the  electrons 

with  spin  a  [30]  and  A  is  the  Anderson  parameter.  Goo,<r(w)  is  the 
local  GF  at  the  QD  site,  which  for  the  immersed  QD  is 

Coo,„M  =  (Gc»)Vcf  M,  (26) 

where  the  localized  Green's  function  Gfa(oj)  is  calculated  employing 
the  atomic  approach,  and  the  conduction  Green's  function  G£(w), 
associated  with  the  leads,  is  given  by  Eq.  (10). 

Following  the  work  of  G.D.  Mahan  and  J.  O.  Sofo  20]  we  define 
the  parameter  f  as  a  function  of  the  thermoelectric  coefficients. 

5  -  ';<T>  (27) 

in  terms  of  what  the  quantity  ZT  can  be  expressed.  Employing  Eqs. 
(22)— (24),  the  dimensionless  thermoelectric  figure  of  merit, 
defined  from  Eq.  (1),  can  be  written  as 


Fig.  3.  The  Wiedemann-Franz  law  WF  vs.  gate  voltage  %D  for  different  hybridizations 
V  and  temperature  T  =  T. 


ZT  =  14f.  (28) 

It  is  clear  from  the  above  equation  that  the  best  dimensionless 
thermoelectric  figure  of  merit  ZT  occurs  at  the  limit  f  — ►  1.  In  the 
next  section,  we  present  the  results  for  the  ZT  and  the  parameter  f . 
We  show  that  in  the  weak  coupling  regime  of  the  SET,  at  temper¬ 
atures  above  the  Kondo  temperature,  it  is  possible  to  obtain  f  values 
close  to  1,  and  in  consequence  the  best  thermoelectric  efficiency. 

5.  Results  and  discussion 

In  our  calculations,  we  employ  the  energy  unit  F  =  1,  where 
D  =  100r  is  the  half-width  of  the  conduction  band  and  the 
chemical  potential  is  always  located  at  /jl  =  0.  The  correlation  U 
generates  a  two-level  system:  one  level  located  at  Eqd  and  the  other 
at  Eqd  +  Uf  which  define  the  behavior  of  the  thermoelectric  prop¬ 
erties  of  the  SET.  The  electronic  correlation  energy  U  is  given  by 
U  =  20.0r,  and  the  symmetric  limit  of  the  Anderson  model  is 
attained  when  Eqd  =  E2  =  -lO.or.  These  parameters  define  two 
other  energies:  E\  =  Ef-  D/2  =  -20.0r  and  E3  =  E/+  D/2  =  0.  Those 
three  energies  are  symmetrical  points  of  the  model.  We  restrict  the 
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calculations  of  the  present  paper,  to  high  temperature  effects.  This 
regime  occurs  approximately  within  the  interval  T=  [l-5]r,  and  is 
characterized  by  strong  charge  fluctuations. 

In  Fig.  2  we  show  the  thermopower  S  vs.  the  gate  voltage  EqD  for 
different  hybridization  values,  at  temperature  T=T.  The  interesting 
point  here  is  the  increase  of  the  thermopower  with  the  decrease  of 
the  hybridization.  The  lower  hybridization  region  (weak  coupling 
limit),  produces  the  most  interesting  results  of  the  present  paper, 
because  the  thermopower  and  the  dimensionless  thermoelectric 
figure  of  merit  ZT  attain  high  values  [29].  This  figure  shows  the 
particle-hole  symmetry  about  the  S_£qd  =  -S£qd  [16].  The  oscilla¬ 
tory  behavior  of  S  is  associated  with  the  competition  between  the 
different  charge  carriers  present  in  the  system:  electrons  and  holes. 
When  S  <  0  the  dominant  contribution  to  the  thermopower  is  due 
to  electrons  and  when  S  >  0  the  dominant  contribution  is  due  to 
holes.  The  figure  shows  strong  variations  of  the  thermopower, 
mainly  in  the  crossover  from  the  empty  dot  to  the  single  occupation 
region  and  in  the  crossover  from  the  single  occupation  to  the 
double  occupation  region.  The  two  oscillations  at  around 
Eqd  =  -lO.or  are  linked  to  heat  transfer  processes  associated  with 
charge  fluctuations.  Recently  an  experimental  paper  [11]  has  been 
published,  that  measures  the  thermopower  of  a  quantum  dot  sys¬ 
tem  in  a  temperature  set  (T  =  0.232  K,T  =  1.5  K,T  =  3.1 1<  and  T=  6 1<) 
closer  to  our  estimated  temperature  (T  =  F  -  10)  I<  for  our  Fig.  2. 
Their  thermopower  results  are  shown  in  their  Fig.  5a,  and  the 
general  behavior  of  the  curves  is  similar  to  that  obtained  by  us  in 
Fig.  2. 

Fig.  3  shows  the  Wiedemann-Franz  law,  in  Lorenz  number  L0 
units  vs.  the  energy  EqD,  at  various  hybridizations  values.  At  low 
temperatures,  the  Wiedemann-Franz  law  is  always  satisfied  (we  do 
not  represent  these  curves  here),  but  at  high  temperatures  it  is  no 
longer  valid,  and  the  departure  from  its  constant  value,  given  by  Eq. 
(20),  shows  a  strong  violation  with  the  decrease  of  the 
hybridization. 

In  Fig.  4,  we  show  the  dimensionless  thermoelectric  figure  of 
merit  ZT  vs.  the  gate  voltage  Eqd  for  different  hybridizations  V,  and 
temperature  T  =  T.  As  the  hybridization  decreases  the  violation  of 
WF  increases  and  for  V  =  4.0r,  it  exhibits  the  highest  departure 
from  its  value  at  low  temperatures.  In  inset  a)  of  the  figure  we  plot 
the  parameter  f,  defined  in  Eq.  (28),  which  exhibits  regions  of 
values  close  to  one,  and  consequently  great  ZT  values.  In  inset  b)  of 


Fig.  5.  Dimensionless  thermoelectric  figure  of  merit  ZT  vs.  gate  voltage  Eqd  for 
V  =  4.0r  and  different  temperatures.  In  the  inset  we  plot  the  parameter  \  vs.  gate 
voltage  Eqd. 


the  figure,  we  plot  the  ZT  values,  WF  and  the  thermopower  S  for 

V  =  4.or. 

It  should  be  stressed  that  the  violation  of  the  WF  does  not  lead 
to  an  automatic  increase  in  ZT  values.  In  Fig.  3,  we  observe  that  the 
maximum  values  of  the  WF  occur  at  around  the  symmetric  point 
EqD  =  -lO.or,  but  the  ZT  values  corresponding  to  this  region,  in 
Fig.  4  are  very  low.  According  to  Eq.  (1),  high  ZT  values  occur 
because  for  some  Eqd  values,  there  is  a  combination  of  lower  values 
of  WF  and  higher  values  of  S  at  the  same  energy  interval  (see  inset 
b)  of  the  figure);  both  properties  contribute  in  a  complementary 
way  to  the  ZT  increase.  The  violation  of  the  WF  relation  is  intimately 
related  to  the  variation  of  the  parameter  f ,  as  indicated  in  inset  b)  of 
the  figure. 

In  Fig.  5,  we  plot  ZT  as  a  function  of  the  gate  voltage  EqD  for 

V  =  4.0r  and  different  temperatures.  For  low  hybridization  and 
high  temperatures  we  obtain  a  strong  ZT  value.  At  the  same  time, 
the  inset  shows  that  in  the  regions  where  ZT  attains  its  maximum 
value,  the  parameter  f  is  close  to  one.  It  should  be  stressed  that  in 
the  ZT  results  presented  in  Figs.  4  and  5  we  do  not  consider  the 
phononic  contribution,  which  tends  to  compete  with  the  electronic 
contribution  to  decrease  the  ZT  values,  as  the  temperature  is 
increased. 

6.  Conclusions 

We  studied  thermoelectric  transport  properties  of  the  single 
electron  transistor.  The  results  obtained  for  the  thermopower  are 
consistent  with  a  recent  experimental  paper  [11].  Quantum  dots, 
due  to  their  sharp  transmittance,  can  be  considered  an  ideal 
nanoscopic  system  for  producing  large  ZT’s  values,  and  one  possible 
technological  application  of  a  single  quantum  dot  may  be  to  create 
localized  cooling  in  nanoscopic  systems  [2].  Quantum  dots  have 
also  been  employed  successfully  for  engineering  quantum  dot 
superlattices.  In  a  recent  paper  Wang  et  al.  [31]  studied  PbSe 
nanocrystal  quantum  dot  superlattices  varying  nanocrystal  sizes 
and  they  showed  that  for  comparable  carrier  concentrations,  PbSe 
nanocrystal  superlattices  exhibit  a  substantial  thermopower 
enhancement  relative  to  bulk  PbSe. 

According  to  Eq.  (21 ),  in  order  to  produce  an  increase  in  ZT,  it  is 
necessary  that  the  increase  of  the  thermopower  S  be  followed  by  a 
decrease  in  the  WF  at  the  same  energy  interval.  We  show  that  this 
condition  is  fulfilled  at  the  weak  coupling  limit  of  the  SET  (low 
hybridizations)  region,  where  the  parameter  f,  introduced  by  G.D. 
Mahan  and  J.  0.  Sofo  [20],  is  close  to  the  unit.  In  this  case  ZT  exhibits 
a  strong  increase.  We  show  that  the  high  values  of  ZT  in  those  re¬ 
gions  are  associated  with  the  violation  of  the  Wiedemann-Franz 
law. 

This  result  shows  that  the  interesting  region  that  could  lead  to 
practical  applications  of  the  system  corresponds  to  temperatures 
within  the  interval  T  =  [1-5]T.  According  to  our  previous  estimate 
in  the  introduction  of  this  paper,  this  energy  interval  approximately 
corresponds  to  T  -  [10-50]  K,  and  can  be  attained  using  a  helium 
liquid  cryostat.  In  this  region,  the  system  shows  strong  charge 
fluctuations,  mainly  in  the  empty  dot  and  the  double  occupation 
regions. 
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